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Using Monte Carlo simulation, we studied fractionalized phases in a model of exciton Bose con- 
densate and found evidence for an emergent photon in a finite size system. The fractionalized phase 
is a meta-stable Coulomb phase where an emergent photon arises as a gapless collective excitation 
of the excitons. We also studied a possibility of spiral* phase where fractionalization and long range 
spiral order of the exciton condensate coexist. 



I. INTRODUCTION 

Recently, possible exotic phases in exciton Bose 
condensate have been studied in the strong coupling 
regime Ij. While the Bose condensation of exciton is 
an interesting subject by itself "2], it turns out that this 
system can serve as a good theoretical laboratory to 
study fractionalization. Fractionalized phase is a many- 
body state where the emergent low energy excitation 
carries fractional quantum number of the microscopic 
degree of freedomQ- Since the Anderson's proposal 
of the resonating valence bond state0, several con- 
crete models have been proposed to show fractionaliza- 
tion, such as exact solvable spin models^, 0], dimer 
model0, bosonic models^, H, 0, 0, and frustrated 
spin 

modelini- 

Fractionalization can be naturally de- 
scribed in terms of deconfinement phase of gauge the- 
ory. The gauge group can be discrete like the Z2 gauge 
theory [a, [a 0] or it can be continuous like the U(l) gauge 
theorvfll IsLlg. ITnl IT^ . While there exist exactly solvable 
microscopic modelsj^, IB| which show the fractionalized 
phase with the Z2 gauge field at low energies, there ex- 
ists no such exactly solvable model or numerical verifi- 
cation for the existence of the fractionalized phase with 
the U{1) gauge field. Therefore it is of great interest 
to numerically show the occurrence of fractionalization 
in the model of exciton Bose condensate which supports 
the U(l) gauge field in the fractionalized phase 1]. Here 
we emphasize that currently there exists no experimental 
realization of the exciton modelQ. We take the model as 
a purely theoretical model to show the existence of the 
fractionalized phase as a proof of principle. 

We consider the Bose condensation of excitons which 
are made of an electron in the a-th conduction band and 
a hole in the b-th valence band where a = l,2,...,Nc 
and b = l,2,...,iVt, denote the band (fiavor) indices. 
In general, the degeneracy of the conduction band Nc 
and the valence band N.^ is different. Here we will con- 
sider the case where Nc — Ny — N. We denote by 
O""^ the phase of the a5-exciton field. When the diago- 
nal (a = b) excitons are condensed, the dynamics of the 
off-diagonal (a 7^ b) exciton becomes relativistic as a re- 
sult of the dynamical constraint posed by the diagonal 



excitons. The off-diagonal exciton is expected to have 
the usual disordered or Bose condensed phases. Interest- 
ingly, we found that in the strong coupling regime a new 
phase can arise where excitons scatter with each other 
to rapidly exchange their constituent particle or holeQ. 
As a result, the exciton loses its identity at a long dis- 
tance scale and a half of exciton which carries only one 
flavor index arises as an effective degree of freedom. The 
emergent fractionalized degree of freedom can be either 
boson or fermion depending on the coupling constants0. 
Although the flavor quantum number of the fractional- 
ized particle is inherited from the original electron/hole, 
the fractionalized particles are different objects from the 
original electron/hole in the sense that they don't carry 
electric charge. Instead they are coupled to a new emer- 
gent gauge flcld. The U(l) gauge theory in (3-|-l)D can 
be in deconflnement phase although the fractionalized 
particles are gapped. The gauge field in the deconfine- 
ment phase behaves in the same way as the photon in our 
world and has been referred to as emergent photon in lit- 
eratures. The gauge field is, on the other hand, nothing 
but a low energy collective excitations of the exciton. In 
the space-time picture, the world sheet of the electric fiux 
line for the emergent gauge field corresponds to the web 
formed by strongly interacting exciton world lines (world 
line web)P|. The boundary of the world line web be- 
comes the world line of the fractionalized particle. This 
explains why the gauge field should emerge as a result of 
fractionalization ! 1 2j . 

It has been pointed out that fractionalization can also 
occur in the presence of conventional long range order 
(dubbed as the AF* phase in the case where the long 
range order is the antiferromagnetic order) because the 
condensation of gauge neutral order parameter field can 
leave the deconfinement phase of the emergent gauge 
theory intact This possibility has also been stud- 
ied in the exciton svstem[T^. Although the previous 
studies^, ^3 provide the rather complete picture, the 
occurrence of fractionalization in exciton system can be 
inferred only in the large TV limit. It is of great interest 
to verify the occurrence of fractionalization by numerical 
simulation at some finite N. This is the purpose of the 
present paper. 
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How do we know if fractionalization occurs ? The most 
direct way would be to probe fractionalized excitation. 
However, it is difficult to directly observe fractionalized 
particles because they are either gapped or strongly in- 
teracting with the gauge field. An alternative way is 
to probe the emergent non-confining gauge field which 
inevitably arises as a consequence of fractionalization. 
The effect of emergent gauge field is most dramatic be- 
cause it is gapless in the fractionalized phase. The ex- 
citon system supports massless U(l) gauge boson in the 
fractionalized phases [HIT^. Currently, there exist other 
models 13, H, IE ^1 which were proposed to support the 
emergent photon in fractionalized phases. However, a nu- 
merical observation of the emergent photon has not yet 
been made in those systems. In the present study, we 
report the observation of the emergent photon in a finite 
size exciton system from Monte Carlo simulation. 



II. MODEL AND PROBE FOR EMERGENT 
PHOTON 

We consider the exciton system in 3-f ID whose Eu- 
clidean action is given by|3, 

^ = E[E{-2«oCOs(0^^^,-0f) 

i a<b 

-2k, E ^os{9'^l^-ef)} 



fj.=2x,2y,2z 



cos I 



qbc 



(1) 



a,b,c 



Here 6f'' = —9\^ is the phase of the off-diagonal exci- 
ton with a ^ b with i, site index in 3-l-lD hypercubic 
lattice. The first two terms are the discrete form of the 
relativistic kinetic energy of the off-diagonal exciton. kq 
and the phase stiffness in the imaginary time and 

spatial directions respectively. In the continuum limit of 

the imaginary time, the ratio — (■^^^ is determined 
from the choice of the discrete time step Ar, where Atq is 
the time step in which the phase stiffness becomes identi- 
cal in space and time. We can study the quantum system 
from the Euclidean action with discrete time step as far 
as we choose the time step such that the temporal corre- 
lation length of exciton is larger than the time step. The 
finite number of lattice in the time direction corresponds 
to the finite temperature of the quantum system and the 
temperature sets the lowest energy scale we can probe 
in the simulation. Taking into account a general time 

step, we set — (■^^) « and k, = (^I') "^here k is 
the intrinsic coupling constant in the isotropic space-time 
lattice. The third term is a coupling between bosons sep- 
arated by two units of lattice in space, which was absent 
in the original model l]. The reason for introducing this 



additional term will be explained shortly. We also de- 
fine the intrinsic coupling constant k, for the frustration 

term as k^ — (■^^^ ■ The last term is the potential 
energy term. Here we consider the cubic interaction with 
> 0. In Ref.,i^ we showed that Eq. ^ is the effective 
action for the exciton model, but it is hard to realize the 
model in real exciton systems for the following reasons. 
First, we are considering the bands with the same de- 
generacy for the valence and conduction bands. Second, 
our result of this paper is that we need either a large 
number of degeneracy or the introduction of frustrated 
interaction for the fractionalization to occur. These con- 
ditions are difficult to realize in real material. Finally, 
we ignore the dipole-dipole interaction between excitons 
which may lead to attraction in three dimension. How- 
ever, the more mathematically minded reader can take 
Eq. as a starting point. It is simply an extension of 
the (d-t-l)D X-Y model to multiple flavors. For exam- 
ple, the effective action may be realizable in Josephson 
junction array. 

In the strong coupling limit = oo, the phase of the 
off-diagonal exciton is constrained to satisfy 

ef = C - <t>l (2) 

where is arbitrary. We will refer to e**^' as slave- 
boson field. Then we can use the slave-boson field as 
our dynamical variables and drop the potential energy 
[K^) term in Eq. Here three remarks are in order. 
First, the replacement of the original exciton field by the 
slave-boson fields is exact in the strong coupling limit 
{K^ = oo) Q . The purpose of introducing the slave-boson 
fields is to parameterize the sub-manifold of finite energy 
in the full phase space of the original exciton fields. Sec- 
ond, the decomposition in Eq. (|2Jl does not necessarily 
imply the occurrence of fractionalization. Fractionaliza- 
tion may or may not occur depending on other coupling 
constants. Third, fractionalization arises only in a long 
distance scale. Microscopically, the slave-boson is always 
confined within an exciton which is the result of infinite 
bare gauge coupling. However, the slave-bosons can be 
regarded as being effectively liberated from exciton in a 
long time (distance) scale as they rapidly exchange their 
partners in the strong coupling limit. We refer to the long 
distance effective degree of freedom as fractionalized bo- 
son in distinction from the microscopic slave-boson. In 
the deconfinement phase, the fractionalized bosons which 
carry only one flavor (band) quantum number arise as low 
energy excitation along with the emergent photon. We 
will search for the emergent photon as a signature for 
the fractionalized phase. The emergence of the mass- 
less gauge boson can be probed from the algebraically 
decaying correlation function of the gauge field. While 
the vector potential of the gauge field is not gauge 

if A dx^ 

invariant, the Wilson operator Wc = e ■'c ^ which 
is defined for a closed loop C is gauge invariant and we 
can express the Wilson operator in terms of the original 
exciton fieldsfi|. We consider the correlation function be- 
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tween the fluctuations of the Wilson operators which is 
defined in the unit plaquctte of the x-y plane, 

Cw{r) = Yl (^W^mU^ + r)SW^2l,i^)) , (3) 

i 

where 

w!:u^) = e^i'^f^'^n.^'^n.^.^'^nj (4) 

corresponds to the Wilson operator of the emergent 
gauge theory and SW^a = W^^- < Kbed >B- V 
is the volume of the system. If the gauge coupling is 
weak the above correlation function measures the flux- 
flux correlation of the emergent gauge theory. In our 
case we expect coupling constant of order unit even in 
the deconfinement phase, and Cw involves higher order 
flux correlation. In any case, an algebraically decay will 
indicate the presence of massless photon and serves as a 
signature of fractionalization. 

III. META-STABLE COULOMB PHASE 

Without K , we found that for iV < 12 a first or- 
der phase transition occurs from disordered to ordered 
phases. These phases are conventional phases where 
there is no emergent photon. Fractionalized phase is 
more likely to occur for a larger N. Here we esti- 
mate a critical Nc for the occurrence of fractionalized 
phase in the absence of the n term based on the mean- 
field approach [J. The reduced action for the slave- 
boson are of quartic form of e'*^ , which can be decom- 
posed by the Hubbard-Stratonovich transformation[l|. 
With the slave-bosons integrated out in the disordered 
phase, the theory reduces to the pure U(l) gauge theory, 
— /3 X]^<iy where is the gauge flux. The inverse 

gauge coupling is given by /3 ~ 2t'^N, where t ~ (N — 
1)kx is the effective phase stiffness for the slave-boson 

with X ' >, the effective hopping inte- 

gral of the slave-boson. Here we used < e''-*» ^ >w< 
eH0?-0-) >< ^-^('P'l-'k)) >^ I < e*W?-<) > assuming 
the flavor independent hopping and ignoring the contri- 
bution of fluctuations. The confinement to deconfine- 
ment phase transition occurs at /3c 1 for the 3-1- ID 
pure U(l) gauge theory^. For k = 0, the order to 
disorder phase transition occurs around tc ~ 0.15 which 
does not strongly depend on N . This implies that the 
minimum number of flavors for deconflnement phase to 
occur near the phase boundary is Nc ^ 1000 which is too 
large to be simulated in a large lattice. In order to reduce 
iVc, we need a larger tc- For this purpose, we introduced 
the frustrated coupling with /t < 0. All of the results 
presented in this paper is for TV = 4. 

With increasing the strength of the frustration, the 
critical tc for the disorder/order phase transition in- 
creases. If the frustration is strong enough, a new or- 
dered phase arises where excitons are 'ferromagnetically' 



correlated within a 2^ box in space and the boxes are 
'antiferromagnetically' ordered with each other. In the 
box ordered phase, the translational symmetry is bro- 
ken. For an intermediate strength of frustration, there 
occurs a box liquid phase where the long range box or- 
der disappears but there still exists strong box correla- 
tion at short distance. This 'box liquid phase' is the 
fractionalized phase where we observe the emergent pho- 
ton. However, it is noted that the region of fractionalized 
phase is very small in the parameter space. To be spe- 
ciflc we set the frustration k N = —0.39 with N = 4. We 
choose a 4^ x 16 lattice with periodic boundary condi- 
tion where the lattice is longer in the temporal direction 
than the three spatial directions. This asymmetric geom- 
etry is chosen because the phase correlation of exciton is 
stronger in the temporal direction owing to the frustra- 
tion in the spatial directions. Further we choose the time 
step so that the effective phase stiffness is smaller in the 
temporal direction than the spatial directions, that is, 
Kx/'^o — (1/0.6)^. With this choice, the strength of near- 
est neighbor exciton correlation becomes comparable for 
time and spatial directions. This effectively enables us to 
study a longer size in imaginary time direction without 
actually increasing the size of lattice. For Monte Carlo 
simulation we used the local Metropolis algorithm. It is 
difficult to implement the more efficient cluster algorithm 
because of the need to introduce the further neighbor 
frustrated interaction. 

In Fig. nia-jj we display the temporal correlation func- 
tion of exciton. For kN < 0.74 the correlation func- 
tion decays exponentially. This implies that exciton has 
nonzero mass gap. For kN > 0.75 there exists a long 
range order and excitons are Bose condensed. The disor- 
dered and the Bose condensed phases are separated by a 
first order phase transition. The spatial correlation func- 
tion also exhibits the first order transition as is shown in 
Fig. [11(b). 

In Fig. 121 we show the mean square of fluctuation of 
the action, < SS'^ >. This quantity is proportional to 
the "specific heat" if we interpret the (3-|-l)D quantum 
system as a 4D classical statistical mechanical system. 
A singularity or discontinuity in the "specific heat" in- 
dicates a quantum phase transition in the correspond- 
ing quantum system. The first order phase transition 
between the disordered phase and the ordered phase is 
marked as a discontinuity in the "specific heat" between 
kN — 0.74 and kN = 0.75. The interesting thing is that 
there is a peak in the "specific heat" between kN — 0.72 
and kN = 0.73, suggesting the presence of another phase 
transition within the disordered phase. In order to un- 
derstand the nature of two different phases within the 
disordered phase we calculate the flux-flux correlation 
function. The temporal correlation function Cw{t) is 
shown in Fig. OKa) and (b). The flux correlation decays 
faster than power-law for kN < 0.72. Since the absolute 
value of the correlation function itself is very small, it 
is hard to determine whether it really decays exponen- 
tially at long distance. However it is natural to expect 
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FIG. 1: (color online) The correlation function between 9^^ , 
the phases of exciton of flavor 1 and 2 in the temporal direc- 
tion (a) and the spatial direction (b). The average is taken 
both for the ensemble and for the site i. The numbers in the 
figure denote the value of kN for each curve. 



that the correlation function of flux operator which is 
composite of gapped exciton fields will also decay expo- 
nentially unless there is special correlation between ex- 
citons. This phase corresponds to the usual disordered 
phase of the off-diagonal exciton where there is no gap- 
less excitations. In the gauge theory picture, this is the 
confinement phase with gapped gauge field. What is in- 
teresting is the algebraical decay of the flux correlation 
in 0.73 < kA^ < 0.74. This algebraic decay of Cw{t) 
is purely due to a nontrivial correlation among excitons 
with different flavors because correlation function of indi- 
vidual exciton is still exponentially decaying. In order to 
demonstrate the importance of the correlation, we com- 
pare the exciton correlation function and the flux corre- 
lation function for kN = 0.73. As is shown in Fig. 0] the 
exciton correlation decays exponentially while the decay 
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FIG. 2: (color online) Mean square of the fluctuation in the 
action as a function of the phase stiffness of the off-diagonal 
exciton. The dashed line identifies a first order transition 
between the disordered phase and the ordered phase. The 
shaded region denote the phase transition within the dis- 
ordered phase from the confinement phase to the Coulomb 
phase. 



of the flux correlation is clearly algebraical. The power- 
law decay of Cw{t) implies that there exists a gapless 
mode. The gapless mode corresponds to the emergent 
photon. This phase is the Coulomb phase where the low 
energy physics is described by the emergent photon and 
the gapped fractionalized bosons Jj. For nN > 0.75, the 
flux-flux correlation function decays exponentially again. 
This is the usual Bose condensed phase of the off-diagonal 
excitons. In this phase, the low energy excitations are the 
Goldstone modes of the excitons and there are {N — 1) 
Goldstone bosons associated with the U{1)^~^ symme- 
try of the model QJ. In the gauge theory picture, this 
phase corresponds to the Higgs phase where the gauge 
fleld becomes gapped due to the Anderson-Higgs mech- 
anism. One of the N slave-boson modes is 'eaten' by 
the longitudianl gauge field and {N — 1) massless modes 
are left. This is consistent with the number of Goldstone 
modes counted from the original exciton model. 

The boundary between the confinement phase and the 
Coulomb phase is not clearly distinguished from CwiT)- 
The flux-flux correlation function decays almost alge- 
braically for kN = 0.71 and 0.72 as well. This can be 
attributed to a cross over behavior at short distance. If 
the conflning scale is larger than the system size, the flux- 
flux correlation function will decay algebraically at short 
distance even in the confinement phase. Although not 
shown here, the absence of a discontinuity in the expec- 
tation value of the action within our error bar suggests 
that the phase transition from the confinement to the 
Coulomb phases is continuous or of a weak first order. 

The observed decaying power of Cw (t) in the Coulomb 
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FIG. 3: (color online) The correlations function Cw between 
the fluctuations of Wilson operators, (a) and (b) are for the 
temporal direction and (c) is for the spatial direction, (a) is 
for 0.69 <kN < 0.74 and (b), for 0.73 < kA^ < 0.76. 
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FIG. 4: (color online) Comparison between the exciton cor- 
relation function and the flux-flux correlation function in the 
Coulomb phase. 



phase is about 2. In the pure U(l) gauge theory, 
the temporal flux-flux correlation function is given by 

< F,,(r,T)^^,,(r,0) >~ Efc,. IJ^Iw^""' ^^^'"^ ^'^■^ 
the field strength tensor and c is the photon velocity. In 
an infinite (3-|-l)D system, the summations become in- 
tegral and the correlation function goes like -p^. The 
discrepancy between the observed decaying power and 
the predicted power may be due to the finite size ef- 
fect associated with the small size in the spatial direc- 
tions. With the small size in the space directions, the 
discreteness of momentum is important and small mo- 
mentum will dominate the momentum summation. In 
an extreme limit, only the smallest possible nonzero mo- 
mentum such as k = (27r/L,0,0) will contribute to the 
correlation function. If the frequency summation is sub- 
stituted with integral, the correlation function behaves as 
— with m ^ 2c7r/L. In reality, the decaying power will 
be between 1 and 4 with an effective 'mass' smaller than 
m ~ 2c7r/L. However, it is hard to make a quantitative 
prediction for the decaying power and the effective 'mass' 
because this analysis is based on the Maxwellian action 
in the continuum. It is expected that non-Maxwellian 
higher order terms of the flux fleld and the contribution 
of large Wilson loops will be also important at short dis- 
tance. Despite the small size in the spatial direction, we 
note that the algebraic decay of the Cw (t) is not due to 
the critical behavior associated with one dimensionality. 
If it were due to the one-dimensional effect, the exciton 
correlation function should show the power-law decay as 
well. 

In Fig. 121(c), the spatial flux-flux correlation function 
is shown. The flux correlation decays more slowly in the 
Coulomb phase compared to the conflnement phase and 
the Higgs phase. However, it is hard to see a power-law 
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FIG. 5: (color online) Mean square of the fluctuation in the 
action as a function of the phase stiffness of the off-diagonal 
exciton. The dashed line identifies a first order transition 
between a disordered (confined) phase and the possible spiral* 
phase. The spiral* phase is a new phase with spiral order in 
9"''' and deconfined emergent boson and photon. The shaded 
region denote the cross-over to the conventional spiral phase 
which is confined. 



behavior of the flux correlation function in the Coulomb 
phase because of the small lattice size in the spatial di- 
rection. The oscillatory behavior of the flux correlation 
in the spatial direction is due to the further neighbor 
frustration k . 
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TABLE I: (color online) The second column represents the 
exciton fields which exhibit spiral long range order in each 
direction and the third column, the slave-boson fields which 
exhibit spiral order in the gauge with (p* = 0. 



IV. POSSIBLE SPIRAL* PHASE 

In view of the finite size effect in the flux-flux correla- 
tion functions, it is desirable to do a systematic finite size 
scaling by increasing the size of lattice and we increase 
the lattice size to 8^ x 16 with N = A fixed. However, in 
the larger lattice it turns out that the Coulomb phase be- 
comes unstable against a formation of spiral long range 
order with pitch 8 in spatial direction. In the smaller 
system with size 4 in the spatial direction, the Coulomb 
phase was stabilized because the spiral order with a pitch 
longer than the system size is suppressed by the periodic 
boundary condition. However, the presence of the long 
range order does not necessarily imply the absence of 



fractionalization as we discussed earlier. In the following, 
we examine the possibility of the coexistence between the 
spiral order and fractionalization. 

In Fig. |S1 we show the "specific heat" with frustration 
K N = —0.3 in the lattice with size 8^ x 16. The time 
step is chosen such that Kx/kq = (1/0.6)^. There is a 
peak in the "specific heat" around kN = 0.54. Although 
now shown here, there is a jump in the expectation value 
of the action between kN = 0.54 and 0.55, signifying a 
first order phase transition. The nature of the phases 
separated by the transition can be understood from the 
exciton correlation function which are displayed in Fig. 
The temporal correlation and the spatial correlation 
in the x-direction show the long range phase coherence 
for kN > 0.55 as is shown in Fig. (a) and (b). On 
the other hand Fig. ini(c) shows that the 1, 2-exciton ex- 
hibits a spiral order in the z-direction for kN > 0.55. The 
pitch of the spiral order is 8. However, it is likely that the 
pitch is not an intrinsic property but is the effect of the 
boundary condition of period 8. For an infinite system, 
it is possible to develop an incommensurate spiral order 
depending on the coupling constants. We proceed with 
the assumption that the nature of the spiral phase does 
not change greatly even though the pitch may change in 
the thermodynamic limit. Similar to the 1, 2-exciton, ab- 
exciton with different a, b develops spiral order in certain 
directions. This pattern of spiral order is summarized in 
Table In the second column of Table we show the 
off-diagonal excitons which exhibit spiral order in each 
direction. It is convenient to understand the pattern of 
the spiral order in terms of the slave-boson field as is 
shown in the third column of Tabled For example, the 
spiral order of 0^^, 9^'^, 9^^ and 6*^^ in the x-direction 
zan be interpreted as the spiral order of (p^ and (fP rela- 
tive to (j>^ and 0^. It is emphasized that the individual 
jlave-boson field can be either phase coherent or inco- 
aerent when the exciton fields have the long range phase 
coherence. In other words, there can be two possible 
phases within the spiral phase. One is the Bose con- 
densed phase of the slave-bosons and the other, the pair 
Bose condensed phase without the Bose condensation of 
the individual slave-boson. The former is the conven- 
tional spiral phase, while the latter is the fractionalized 
phase with emergent photon. The latter possibility has 
been discussed in the uniform condensate called Higgs* 
phase^^. It can be understood as the condensation of 
bundles of vortices of individual slave-bosons without the 
condensation of vortices of the individual slave-bosonfl^. 
Similar consideration applies to the spiral phase and in 
order to distinguish the conventional spiral phase from 
the fractionalized spiral phase, we will refer the latter 
phase as spiral* phase. 

Unlike the Coulomb phase, it is very difficult to unam- 
biguously identify the spiral* phase, because other gap- 
less fluctuations are present in the form of iV — 1 Gold- 
stone modes. It is necessary to minimize the contribu- 
tion of the Goldstone modes to the flux-flux correlation 
function. For the reason we calculated the flux-flux cor- 
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FIG. 6: (color online) The correlation functions between 9^^ , 
the phases of exciton of flavor 1 and 2 in the temporal di- 
rection (a), the a;-direction (b) and the z-direction (c). The 
average is taken both for the ensemble and for the site i. The 
numbers in the figure denote the value of kN for each curve. 



relation function defined as 

C'wir) = -Re^ E {^Wi2ti^ + r)5W^U^)) , (5) 

i 

where VF^j'^^(z) is the Wilson loop operator defined in 
Eq. Q). Here we consider the flux- flux correlation func- 
tion which is different from Eq. and the correlation 
is calculated between two VF^^"^^ whose flavor indices are 
ordered in a reverse wayj^fij. This is to suppress the con- 
tribution of Goldstone modes in the presence of the spiral 
long range order. The Goldstone mode comes from the 
U global symmetry under which the phase of exci- 
ton transforms as 0°'' — > Of'' + (p"' — cp''. Since the Wilson 
operator does not carry a net flavor, it is coupled to the 
Goldstone modes only at a flnite energy /momentum as 

However, the first derivative term in the above expres- 
sion does not contribute to C^^(r) because of the differ- 
ent ordering of flavors in the two Wilson operators inside 
C^y(r). In energy-momentum space, the first derivative 
terms contribute to the temporal correlation function as 
/ d^kk^ky < ipl(p% > e^'^-o^ - / d'^kk^ky/k'^e^^°'' and it 
vanishes. As a result, the first non- vanishing contribution 
from Goldstone modes come from the second derivative 
term and C'^{k) is coupled to the Goldstone mode with 
an extra factor of fc** (In contrast the correlation Cw in 
Eq. ||2J) will only have a factor of fc^) . Upon Fourier trans- 
formation, the Goldstone mode contributes to C'^{t) as 
l/r^. If we take into account finite size effect of spa- 
tial directions, C^y(r) can decay at most as Thus 
the power-law decay of the temporal correlation func- 
tion with a power smaller than 3 can not be explained in 
terms of Goldstone mode. On the other hand, the emer- 
gent photon mode will contribute to C^y(r) as in 
an infinite lattice just like the fiux-flux correlation func- 
tion in 3+ ID electrodynamics. The flnite size effect can 
modify the behavior upto 1/r. 

In the spiral phase the temporal flux-flux correlation 
function decays algebraically as is shown in Fig. 13 (a) and 
(b) . Result in the disordered phase is not shown because 
it is too small (less than 10~^ at nearest neighbor) and 
further neighbor correlation is completely dominated by 
statistical noises. The decaying power of Cy^{T) is about 
1.3 for kN = 0.55 and it increases with increasing k. 
The decaying power becomes about 3.3 at kN = 0.9 
and 3.6 at kN =1. It first becomes greater than 3 
around kN = 0.85. Thus the spiral* phase may exist 
for 0.55 < kN < 0.80 and the conventional spiral phase, 
for kN > 0.85, with a continuous phase transition in be- 
tween. However, it is a possibility that the slower decay 
for kN < 0.80 is due to a cross-over behavior within the 
conventional spiral phase. This is because the "specific 
heat" does not exhibit an anomaly within our error bar 
for the phase transition from the spiral to spiral* phases 
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as is shown in Fig. |S1 In the spatial flux-flux correlation 
function, the algebraic decay is less clear because of the 
finite size effect. However, it is clear that the correlation 
decays more slowly for smaller k as is shown in Fig. [3(c). 
Although not shown here, the flux-flux correlation func- 
tion has the same dipolar properties as < hz{r)hz{0) >, 
where bz is the z-component of the 'magnetic' field in 
the emergent gauge theory. As expected from the gauge 
theory, Ct^ (r) is positive (negative) if r is parallel to the 
z(x) -direction. 



V. CONCLUSION 

In conclusion, we observed the emergent photon in a 
meta-stable fractionalized phase of the exciton Bose con- 
densate. In the Coulomb phase, the photon mode arises 
as a gapless collective excitation although individual exci- 
tons are gapped. The Coulomb phase is only meta-stable 
because it is unstable in a large lattice and a spiral phase 
becomes stable. In the spiral phase we studied a possibil- 
ity where fractionalization coexists with the long range 
spiral order. We observed a slow collective mode of ex- 
citon which can not be explained in terms of pure Gold- 
stone modes in the spiral phase and interpret the slow 
mode as the emergent photon in the fractionalized spiral 
(spiral*) phase. 
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FIG. 7: (color online) The correlations function C^^ between 
the fluctuations of Wilson operators, (a) and (b) are for the 
temporal direction and (c) is for the z-direction. In (a), dif- 
ferent scale factors are multiplied to separate the curves as 
denoted in the figure. In (b), the correlation functions are 
shown for k,N — 0.55 and 0.9 with bare scale. 
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